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ABSTRACT 


In this report, the analysis, design and operation of single-phase, single-slot tubu¬ 
lar permanent magnet linear alternator is presented. The study includes the no-load 
and on-load magnetic field investigation, permanent magnet’s leakage field analysis, 
parameter identification, design guidelines and an optimal design of a permanent 
magnet linear alternator. For analysis of the magnetic field, a simplified magnetic 
circuit is utilized. The analysis accounts for saturation, leakage and armature reac¬ 
tion. 




List of Principal Symbols 

B Flux density, T 

f m Frequency of the alternator, Hz 

H Magnetic field intensity, A/m 

/ Demagnetizing component (of armature reaction) of stator (or primary) 
current, A 

/, Stator (or primary) current, A 

k A coefficient 

k M Fill factor 

k gm Leakage coefficient of permanent magnet. 

L„ Leakage inductance, H 

L m Mutual inductance, H 

L me Equivalent mutual inductance, H 

R„ Stator (or primary) winding resistance, n 

R l Load resistance, n 

W Weight, lb. 

Wi Stator (or primary) winding turns, turn 
Z { Iron loss, watt 

Z c Copper loss, watt 

<j> Flux, Wb. 

X Flux linkage, Wb-turn 

w Electrical angular frequency, R/S 

Subscripts 

gav Average value in airgap 

pch Average value in the horizontal part of flux path in plunger back iron. 

pev Average value in the vertical part of flux path in plunger back iron. 

sp Average value in stator pole shoes. 

tcv Average value in the vertical part of stator back iron. 

sch Average value in the horizontal part of stator back iron. 

g airgap 

p plunger 

9 stator 

m magnet 

v vertical 

h horizontal 

c core (back iron). 
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CHAPTER 1 


INTRODUCTION 

Because it is a self-starting, quiet vibration-free and highly reliable device, the 
free-piston Stirling engine (FPSE) has recently received considerable attention. The 
FPSE is a strong candidate as a prime mover to drive a linear alternator in space 
power stations. The FPSE may use solar energy or nuclear energy as the source and 
can perform in any orientation—vertical, horizontal, inclined or upside down. For 
coupling with the FPSE, the linear alternator (LA) is most suited because of its sim¬ 
ple structure and high efficiency of direct coupling. Fig. 1-1 shows a 5 kW free pis¬ 
ton Stirling engine—linear alternator unit. 

According to the method of excitation the linear alternators may be: (i) field- 
excited, with moving diodes; (ii) inductor or reluctance type, having both armature 
and field windings on the primary or (iii) permanent magnet (PM) excited. The PM 
linear alternator does not need any exciting winding and exciting power supply, so 
that it is simpler more reliable and lighter than the field-excited or reluctance-type 
linear alternator. Hence the PM linear alternator is considered more suitable for 
space power applications. 

Past work on PM linear alternator is as follows: 

(i) Primitive PM linear alternators are presented in [l], [2] and [3]. 

(ii) References [4] and [5] discuss the fundamental equations and basic design guide¬ 
lines of a PM linear alternator. 

(iii) Reference [6] gives one design of PM linear alternator, but the design is not 
optimal for either the ratio of power to weight of the alternator or the efficiency 
of the alternator. 

It is fair to say that the above does not reflect adequate analysis and design of 
PM linear alternators. The present work consists of the field analysis, parameter 
identification, design guideline and optimal design of PM linear alternators. 

The following discussion is for the single-phase, single-slot tubular linear alterna¬ 
tor. But the results can be extended to multi-phase and multi-slot tubular linear 
alternators. 
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CHAPTER 2 

THE FIELD OF A TUBULAR PM LINEAR ALTERNATOR 


2.1. Configuration of Single-Phase Single-Slot Tubular PM Linear Alter¬ 
nator 

A single-phase single-slot tubular PM linear alternator is shown in Fig. 2-1, 
which has four pieces of magnets mounted on the plunger. The circular stator coils 
are embedded in stator slot. The plunger reciprocates with the stroke length: 

t*roke = (2-1) 

The flux 4 > of permanent magnets in the airgap varies almost linearly with the 
position of the plunger, as shown in Fig. 2-2. If the plunger moves with a sinusoidal 
speed, a sinusoidal voltage will be induced in the stator winding. The flux of the per¬ 
manent magnets has a maximum value at the ends of the plunger stroke. 

2.2. Basic Equations of Magnetic Field and Demagnetization Curve of 
Permanent Magnet 

Two very useful laws for electromagnetic field analysis are Ampere’s Law and 
continuity of flux. Let there be n parts in a magnetic loop, then by the Ampere’s 
Law: 


WI= £ HA (2-2) 

i-i 

Continuity of flux requires that 

= <t>2 = ' ’' (2-3) 

In (2-2) and (2-3) W is the coil turns linking the magnetic loop, I is the current in the 
coil, fa, l { are field intensity, flux and length of part i of the magnetic loop respec¬ 
tively, i=l,2,...n. 

Let TASCORE 27 be the permanent magnet material used in the linear alterna¬ 
tor for which the demagnetization curve is shown in Fig. 2-3. It may be seen that 
the demagnetization curve of TASCORE 27 is almost a linear curve and may be 
described approximately by: 


H m = H e ( 1 - ^-) 

(2-4) 

B m = B r (1 - -^) 

(2-5) 





Figure 2-1 Linear alternator geometry. 
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Figure 2-2 Flux variation with plunger position. 
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2.3. Procedure to Calculate No Load Field 


The dimensions shown in Fig. 2-1 have the following relationships: 


A* - A -2 9 -h m 

(2-6) 

A-1 = A -2 g - 2 h m - b p 

(2-7) 

D P 2 - A -2g-2h m -± b p 

(2-8) 

A = A 4 hy 

(2-9) 

Ai = A* + + A 2 + y 

(2-10) 

A 2 “ A 4 2Aj 4* 2A 2 4 6* 

(2-11) 

b = 6# 2 4- frf 

(2-12) 

The magnetic loop in Fig. 2-1 is divided into several sections, 
that the field is constant in each section, and by applying Ampere’s 
ABCDEFA, we obtain: 

If it is assumed 
Law to the path 

+ Ac*(2A2+M 4 Ap^Ai 4 


4 + H m 2h m = 0 

(2-13) 

or 


Hm= ~2h^ K * 6 + ^’ (2A2 + h) + H,p2hl + Hg ” 29 


+ H pep b p + H pek ij 

(2-14) 

It will be shown in the next chapter that the leakage flux of permanent magnet 
is not a constant and can be expressed as: 

4>, = -kH m =-kH e (l-^) 

(2-15) 

where k is a constant determined by the dimensions of the PM alternator. Con¬ 
tinuity of flux (Fig. 2-4) yields: 

r m D m B m + kH e (1 - £ ) - + 2 m (A - 9) B g „ 

(2-16) 

L Y !!L (A - 9) Ar« - 7 (r + b ( ) D,B, P 

(2-17) 

“ (r 4 &*) AAp = b t D n Ac* 

(2-18) 

b t D tl B $cv = ifAsAcA 

(2-19) 

An An — D p2 B pc9 

(2-20) 

r mB p nB pcv = b p D pi B pch 

(2-21) 

The saturation curve of the core material (M-19) is shown in Fig. 2-5. The 


algorithm for the determination of the field of the permanent magnet is as follows: 
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Figure 2-4 PM flux paths and continuity of flux. 
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(i) Choose a value of H mo on the demagnetization curve and an accuracy criterion t. 

(ii) Find B mt from (2-5). 

(iii) Find B in every part of the magnetic circuit with B mo and (2-16) through (2-21). 

(iv) Determine H for each B from the B-H curve (Fig. 2-5). 

(v) Determine H g „ = B g „/n 0 for the airgap. 

(vi) Calculate H ml from (2-14). 

(vii) If 1 / 4,0 - H ml I < stop calculations. The field due to this H m0 is the field of the 
permanent magnet. 

(viii)If I H m0 - H ml I > €, choose (H mo + H ml ) as a new H mo and return to stop (ii). 

2.4. Field of a Linear Alternator on Load 

The flux distribution in a linear alternator is shown in Fig. 2-6, where <f> x denotes 
the mutual flux due to the permanent magnet, the mutual flux due to the stator 
mmf, <j> = - fa the net mutual flux and 4>„ the leakage flux due to permanent mag¬ 

net. For ease of calculation, certain cross-sections are defined as follows. These 
cross-sections are denoted by S with subscripts, 


Sgl = * (A -»)“(»■ + T m) 

(2-22) 

s t u - x (A - 9){t + 2 [g + 9m)—^ 7—] - *(A - 9){t + 9 ) 

9 + n m 

(2-23) 

Sg2 = X (A - 9 - A m)( r + 9 + Aj 

(2-24) 

Sgi = x ( A -2 g - A m )(r + 2 g + h m ) 

(2-25) 

Sg3m - X ( A -2 g - h„) T m 

(2-26) 

5 v -t(A + A 1 ) j(r+b t ) 

(2-27) 

S Kt = jt (D t + 2 A] + A 2 + — b k ) b t 

(2-28) 

S*ck 38 x (A + 2 A t + 2 A 2 4- it) b b 

(2-29) 

Spc, - X (A - 2$ - 2 A m - J 6 p )(r + 2 ? + 2 A m ) 

(2-30) 

*^pcA “ T (A ~ — 2 A m — 4 p )4 p 

(2-31) 


Let 4> m denote the flux due to the permanent magnet within the permanent mag¬ 
nets. Then: 


<t>m = <t>a + ^1 

(2-32) 

and in the permanent magnet: 


u 

1 3 

1 

(2-33) 

In the airgap: 


o „ (? + *») - - AnAm (— —) 

(2-34) 





Figure 2-6 Fluxes in a LA 
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In the iron portion of stator, total flux 


11 

1 

& 

(2-35) 

ii 

(2-36) 

B = 

(2-37) 

D _ <t> 

t >$ch - 

^•ck 

(2-38) 

In the iron portion of the plunger: 


n o 

^pct — c 

(2-39) 

°pci 


r, 

Vpch - ~ 

Jpek 

(2-40) 


If <j> u 4> 2 and <t>, are known, flux densities in every portion of the iron can be 
obtained from (2-32) and (2-35)-(2-40). For the corresponding flux density, the field 
intensity can be found from the magnetization curve of the core material. 

Let: 


F m = 2 H„h x + 2 H,„ (A 2 + I M + H Kk b + 2 H pe , (| b p ) + H pck b 

(2-41) 

By Ampere’s Law: 


(f)r> 20o 

, ‘ F w + ; (? + A m ) - wy 

(2-42) 

From Eq. (2-4), (2-15), (2-32) and (2-33): 


rr rr H ‘ f^l ~ kH m ^2 , 

n L c C J 

d t dgZ 

(2-43) 

or, 


i, m - in ^ n 4 

[1 DC 1 D C D C 

&r^g3 m &r^qZ 

(2-44) 

Let 


„ kH e 

d = i -‘-r- 
BrSg3m 

(2-45) 

Equation (2-44) becomes: 


,rr rr H ^l , 

dM m *i c dc d c 

(2-46) 

and 


*1 = [(** - dHJ 

(2-47) 

Substituting (2-47) into (2-34) yields: 
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Or, 


Let: 


and 


Now 


---(* + U = -H m h m - l(H e - dH m ) f- + f-] 

s g2^o M c ^gZ ^glPo ^gll^o 


B r Sg3m9 

Ni 


PtK ~ 


g + N3mg _ g 

S t 2 SglSgZ Sg\\ 


^BfSgZm g 

Ml 


g + Ni SgimQ 

Sg2 Sgl$g3 $gl 1 




BrSgim 9 


a — 


_£L 


9+_Kn Sg3m9 0 

rt ^ n t-t n 


J D2 


NlN* 


/'oN, - 




dBfSgim g 

Ml 


g + ^ ^3mg _ g 


^2 


Nl^l/3 Nil 


02 — « + 

Substituting (2-52) into (2-42) gives: 

0 2 Fm + *L±M {a + ^ = ^ 


01 — 02 
- 1 (- 


Sg2^0 




(2-48) 


(2-49) 


(2-50) 


(2-51) 


(2-52) 

(2-53) 

(2-54) 


To find the field of the linear alternator, choose (or guess) an initial value H mo 
(of H m ). From (2-52), (2-47), (2-15), (2-32) and (2-35)-(2-40), the 0 2 , 0 1? <j>„ and the flux 
densities in all parts of iron corresponding to H mo can be found. Another value of H m 
is obtained by the magnetization curve of iron and Eq. (2-41) and (2-54). After a few 
iterations exact values of H m and the fluxes are obtained. 


2.5. An Example 

A single-slot single-phase tubular PM linear alternator has the configuration 
shown in Fig. 2-1 and has the following dimensions and parameters: 

r = 1.2"; r m = 1.0"; A, = 0.2"; h 2 = 0.8"; g = 0.036"; 

h m = 0.46"; b t = 0.5"; b b = 0.45"; b p = 0.7"; b, 2 = 3.44"; 

D { = 11"; B r = 1.07 T) H e = -770.305 kA/m; W 1 = 21 turns 

By the calculating algorithm in section 2.3, the no load field of the alternator is 
calculated. The B and H in every part of the machine are: 

B m = 0.97242 T H m - -70251.93000 A/m 
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B, p = 0.93418 T = 81.77248 A/m 

B, ev = 1.43153 T H k , = 973.84080 A/m 

= 1.46938 T H uh = 1698.5440 A/m 

£ pe , — 1.05397 T H pev = 100.79380 A/m 

B pck = 1.56229 T = 3477.81500 A/m 

B ga , = 0.73740 T //,„ = 585241.40000 A/m 

The mutual flux produced by PM is: 

<t> = 0.01817 IW 

The leakage coefficient of PM is: 

Km = 1 + - 7 - - 11446 
9 

When the stator current /, = 209.3 A and 0 = 45* (shown in Fig. 2-7). The 
demagnetizing component of the stator current is: 

I = 7[sin0 = 209.3 X sin 45" = 148 A 
By the formulas given in section 2.4, the field is: 

B m = 0.85001 T H m = -158371.40000 A/m 
= 0.50827 T H, p = 40.57893 A/m 
B,„ = 0.77887 T H, ( , = 63.09837 A/m 

B, ck = 0.79946 T H Kh = 64.95149 A/m 

B pe , = 0.80115 T //p C , = 65.14345 A/m 

5 p a = 1.18753 T H pch = 127.50650 A/m 

The mutual flux due to PM is: 

4> x = 0.01549 Wb 

The mutual flux due to the current I is: 

<t> 2 = 0.00561 Wb 

Net mutual flux: 

0 = < 0 , - <t> 2 = 0.00988 Wb 

The leakage coefficient of PM is: 

K m = 1 + — - 1.3822 
a \ 

The coefficient k am on load is larger than that for no load. The reason is that 
the magnitude of H m increases with an increase in I. 
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CHAPTER 3 


THE LEAKAGE FIELD OF PERMANENT MAGNETS 
IN A PM LINEAR ALTERNATOR 


The leakage of the permanent magnets in a tubular PM linear alternator is 
shown in Fig. 3-1. Here the symmetrical line of a permanent magnet coincides with 
that of the stator teeth. For other positions of the permanent magnets the leakage 
field may be found in a way similar to the following. The field shown in Fig. 3-1 is 
only in the vicinity of a stator pole. The entire region of the field is divided into a 
number of sub-regions as shown. Flux calculations are made for these sub-regions. 

In the flux calculations, the following conditions are assumed: 

(a) In the sub-region ABCDA, the fluxes are perpendicular to AB and CD. 

(b) In the sub-region BCGFEB, the flux pathes are semi-circles with different diam¬ 
eters and the center of those circles is the middle point between B and C. 

(c) In the sub-region HIJKH, the flux pathes are quarter circles with center at point 
H. 

(d) In the sub-region IJKPOMLI, the flux pathes are composed by quarter circles 
with center at point H and quarter circles with center at point K. 

With the above assumption, in sub-region ABCD: 


F m = -H m h 


dK 


_ Frn _ 

_ gJ2 _ 

HoA D i -2 g - 2 h m + 2h)dh 
-2-Hm^ a ?r(A- ~2g - 2A m + 2A)A ^ 
9m 


(3-1) 


(3-2) 


d><n = / 



-2-2 g_- 2h m + 2ft 
9m 


dh 


9m 


[(A -2? 


2A m ) A* 



In sub-region BCGFEB: 

r, - A/r + (j (Z + 7 , - ~ 2 

F.-- 2H m h„ 

The reluctance on a flux path with width dr is: 

in _ _ T_ 

2 - g)ii„dr /»,(A ~ 9)dr 


(3-3) 


(3-4) 

(3-5) 


(3-6) 
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Figure 3-1 Leakage field of the permanent magnets. 
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where A - g is the average diameter of the cross-section of the flux path. 
The infinitesimal flux over the path is: 

= = - 2 ^H m h m (A -?) — 

an 2 r 

r i 

$o2 “ / <^«t2 

tfm/ 2 

= —2H m fi 0 h m (D; - ?)ln — 

9m 


In the sub-region HIJKH: 

A. = -A. A 

jrA_ 

<f/?3 = tt(A - 2g - A m ) = 2(A - 2? - A m ) J 0 Ih 
where A - 2? - A m is the average diameter of cross-section of the flux path. 

^03 = - - ^An^e (A ~ 2 ? - A m ) (to 

or 

^<73 — / <^<73 
0 


(3-7) 


(3-8) 


(3-9) 

(3-10) 


(3-11) 


= -2H m fx 0 (D { -2g - A m ) h m 

In the sub-region of LMOPKJIL: 

F m = -H m h m 

±*(h m +h) + ±xh 

dR i = ^-- - 7 r~~7 

*(A ~ ~ A m )/i„(/A 


A m 4- 2A 


2(A - 9 ~ A m K rfA 


^ ~ rfi? 4 


Pm —2A m /i„(A - 9 - A m )A„ 


A m + 2A 


rfA 


^<74 — J ^4>CTi 

0 

= -A, (A - ? - A m )ln ) A n 


Total leakage flux: 

^17 = 0<7l + ^<7l + 0(73 + 0(74 


(3-12) 


(3-13) 

(3-14) 


(3-15) 
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HmPo 


7T 

9m 


A - 2 ? - 2A m )A m + |- A* 


+ 2(A - ?)ln 




+ 2(A ~ 2? - A„) + (A- - 9 - A„)n( 


,*m + 2y, 


where r t is given in (3-4). 

The Eq. (3-16) may be rewritten as: 

= -KH m 

where 

+ 2(A - ?)ln —- 

9m 




(A - 2? - 2A m )A m + -J A* 


+ 2(A - 2(7 - A m ) + (A - (7 - A m ) In 



A m 


(3-16) 


(3-17) 


(3-18) 


It is obvious that K is dependent on the configuration of the alternator. The 
flux <b m in the permanent magnet is given by: 


K = * (A - 2(7 - A m ) r m B m 


(3-19) 


Combining (3-17), (3-19) and (2-4) yields: 
4> a = -K(H C - jj- B m ) 

H c i 


= ~K H, + K -rr~ 


B r ?r(A -2g - A m ) r„ 


(3-20) 


Equation (3-20) shows that the leakage flux decreases linearly with B m (or This 
happens because 4>„ is directly proportional to H m , which decreases linearly with B m . 
The conclusion is that after the dimensions of the linear alternator have been 
decided, to reduce leakage flux B m has to be chosen as high as possible. 

The above analysis may be used to calculate the leakage coefficient of the linear 
alternator. By definition the leakage coefficient is: 


k 


am 


K 


(3-21) 


From (3-20): 


k = 


K +KH t -K 


H, 


B r jt(A - 2(7 - A m )r m 


1 


1 +KHA-J-- 


_i_ 

B r n {D { -2 g 



(3-22) 
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Defining: 


and substituting (3-19), (3-23) into (3-22) yields: 


Km - i -r- 

1 + kH c ( 1 - k B ) -f- 

rm 


(3-23) 

(3-24) 


Since H c is negative, when B m approaches B r , i.e. k B approaches 1, k am decreases 
to the limit 1. If B m decreases, because 4> m will decrease with B m , the second term in 
the denominator of (3-24) will increase rapidly. Thus k am will also increase rapidly. 
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CHAPTER 4 

PARAMETERS OF PM LINEAR ALTERNATOR 

The three important electrical parameters of a single-phase PM linear alternator 
are: leakage inductance, magnetizing inductance and resistance of the stator wind¬ 
ing. The methods and formulas to find the three parameters are presented as fol¬ 
lows. 


4.1. Determination of Slot Leakage Inductance 

The shape and the dimensions of the slot of a tubular single-phase single-slot 
tubular PM linear alternator are shown in Fig. 4-1. The winding in the slot has W x 
turns. The slot is divided into two regions—ABCF and CDEF as shown in Fig. 4-1. 
Let the current in the winding be 

For calculating the slot leakage inductance, the following assumptions are made: 

a) The current density in the coil is uniform. 

b) Iron around the slot has a much larger permeability than that of air. 

c) The leakage flux in the slot is parallel to the bottom of the slot. 

In region ABCF, over a depth dx, the infinitesimal flux linkage is: 

d \, = ( W, ff *(D, - 2x) dx (4-1) 

"2 °»2 


The total flux linkage in the region ABCF becomes: 


Xi — J* — 


P,*IxW? 


6,2 i v h 


J(-*-f(D-2x)dx 




(4-2) 


For the region CDEF, the approximate equivalent slot width is — (6„ + 6, 2 ). 


Thus: 


</X 2 — wh w? (■ 


6,1 + 6, 2 


') (Ai + 2 y) dy 


X 2 = / d\ 2 = fi 0 nIiWi 
o 




6«i + 6, 2 0 


(A*i + h-) 


/ (A + 2 y) dy 


6,i + 6, 2 

Total slot leakage flux linkage is: 

X,„ = X t + Xj, 

Hence, from (4-2), (4-4) and (4-5), the slot leakage inductance is: 

, _ ^ _ lv2 , A fh 2 AA 2Af 

,a I x 1_ 36, 2 26,2 6,i + 6,2 6,i + 6, 2 


(4-3) 


(4-4) 

(4-5) 

(4-6) 
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where 


D 0 — D{ + 2A] + 2Aj (4-/) 

Similar expressions may be obtained for other slot geometries. 

For a traditional rotating alternator, with same shape of slots as shown in Fig. 
4-1, the leakage inductance of one slot is: 





2A' 
b>i + 



(4-8) 


where l ef is the length of slot. 

Comparing Eq. (4-6) and (4-8), there are two additional terms—term 2 and 
term 4—in (4-6). The reason is that unlike a traditional rotating alternator, the slot 
length of a tubular alternator is not the same for different diameters in the slot. 


4.2. Determination of the Magnetizing Inductance 

Applying Ampere’s Law to the path ABCDEFA in Fig. 4-2, with a primary coil 
current I u we obtain: 

H. ch (b) + H kv { 2 A 2 + b h ) + A/^( 2 Aj) + H 9 „2(g + M 

+ HUh) + H Pck(b) = F « W X I X (4-9) 

In Fig. 4-2, the flux distribution is given. With the dimensions defined in (2-6)- 
(2-12) and applying the continuity of flux condition across the five set of cross- 
sections in Fig. 4-2, we have: 


( r + 9 + A m )(r> m + g) B g „ = -J- (r + b t ) D,B, P 

(4-10) 

~ ( r + b t )D, Bfj, = b t D n B, ev 

(4-11) 

btDxxBw — b t D t 2B Kh 

(4-12) 

(r 4- h m )(D m -h g) B^ a9 “ (r -f 2g 4- 2A m ) D p2 B pc9 

(4-13) 

(r + 2g 4- 2h m )D p2 B pcv = b p D pl B peh 

(4-14) 


For a certain value of B t „, the flux density in other parts of the magnetic loop 
of the linear alternator can be obtained from Eq. (4-10) through (4-14). The 
corresponding field intensities are obtained from the magnetization curve of the core 
material. For air: 

H g „ = j- B t „ (4-15) 

After the field intensities in every part of the magnetic loop have been found, 
the magnetomotive force F needed to maintain the value of B ga , could be calculated 
by Eq. (4-9). The flux linkage, X m , is given by: 

x m = W x <t> = W x B g „x ( D m + g)(r + g + AJ (4-16) 

Hence, the magnetizing inductance, L m , may be expressed as: 

.... W?£ | ;i» ;r i£m + g)( T + 9 + ^m) 

Uvh) — I — r. 


(4-17) 
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C-C'r Q, 



Figure 4-2 Flux path and continuity of flux. 





- 25- 


Notice that this inductance is a function of the primary current Thus, the effect 
of core saturation is taken into account globally and L m depends on the value of B g „. 
Fig. 4-3 shows the variation of the L m with B g „ in a tubular linear alternator. 

4.3. The Equivalent Magnetizing Inductance of PM Linear Alternator 

The method discussed in the preceding section is useful in finding the magnetiz¬ 
ing inductance of a tubular linear alternator, such as self-excited, reluctance as well 
as the PM type. But for a PM linear alternator, the magnet’s working point changes 
with the variation of load. Thus, the induced voltage will be reduced not only by the 
current reaction but also by the variation of the working point of the magnet when 
the load increases. For design and analysis of PM alternators, sometime it is con¬ 
venient to define an equivalent magnetizing inductance as follows: 

«- 18 > 

where 

E 0 = no load induced voltage 

E = airgap induced voltage corresponding to stator current I x and angle 0 shown 
in Fig. 4-4. 

w = electrical angular frequency. 

Now, the equation of circuit of Fig. 4-4 is: 

E„ — jwL me /i + ]uL , 9 fj + R t I { + V (4-19) 

Consequently, L mt {I x ) is called an equivalent magnetizing inductance. 

4.4. Resistance of Stator Winding 

The resistance of stator winding is given by: 

Ra - + 2*i + An) (4-20) 


where 

p = the resistivity of stator winding material 
A e - cross-section of stator winding wires 
D { + 2A t + h 2 = average diameter of stator winding 
IF, = turns of stator winding 

4.5 An Example 


For the PM tubular linear alternator, with the dimensions given in section 2.5, 
the leakage inductance of the slot is: 




(Z?,’ + 2Aj + 2 Ao)A 2 Ao 


2 A A. 


2Af 


3 6 


<2 


2 i,2 &,1 + A,2 6,1 + 6,2 


(11 + 2 X 0.2 + 2 X 0.8) X 0.S 


0 . 8 2 


3 X 3.44 


2 X 3.44 


= 7T X 21 2 X 1.26 X 10-* 



< MW > 
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B (T) 


Figure 4-3 The magnetizing inductance versus flux density 
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2 X 11 X 0-2 2 X 0.2 2 

2.04 + 3.44 + 2.04 + 3.44 


X 0.0255 


= 0.0771 mH 


For /j = 210 A and 9 = 45°: 

Using the method introduced in section 4.2, the magnetizing inductance is: 

L m = 0.822 mH 

By the definition given in Section 4.3, the equivalent magnetizing inductance is: 
L mt = 0.832 mH 

The constant k, which is calculated by (3-18), is: 
k = 0.3738 X 10~ 7 Wb m/A turn 

From this example, it can be seen that the value of L m is close to that of L me . 
But if k is much larger than the above value for designed alternator, the difference 
between L m and L mc will be large. 
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CHAPTER 5 

DESIGN OF A TUBULAR PM LINEAR ALTERNATOR 

In Reference 5, the basic design guide of PM linear alternators with tuning 
capacitor is given. The design procedure of PM linear alternators without tuning 
capacitor is presented as follows with a design example. The configuration and 
dimension symbols of the designed alternator are the same as those shown in Fig. 2- 
1. 


5.1. Given Data 

Rated output power (s) = 25 kVA 

Frequency (/ m ) = 100 Hz 

Plunger travel length =1.0 inch 


Assumptions and Estimates 

Length of pole shoe of stator r = 1.2 in 
Pole length of plunger, r m = l.o in 
Airgap g = 0.036 in 

Full-load output voltage V = 120 V (rms) 

Induced (or internal) voltage E - V2(120) V (rms) 

Maximum induced voltage E m = 2(120) = 240 V (max) 

Output current 7=4= 25 ^ = 208.33 A 
* V 120 

Magnet material = Tascore 27 H B r = 1.07 T, H e = 770.305 KA/m 


5.3. Design Procedure 

Let B m = 0.97 T (satisfied after iterations) 


77 m = H c {1 — ■“*) = 


B ga , = 


B m r m 


~ (r m + r)K a 


770.305 (1 - 

_ 0.97 X 1 
" 1.1 X 1.05 


0 97 

= - 71.99 kA/m 
1.07 ’ ' 

= 0.7688 T 


where k am = leakage coefficient = 1.05 (assumed) 
H m h m = ^-gk, = H t „ g k, 

r^o 

where k, = no-load saturation factor = 1.5 (assumed) 
From (5-3): 

A _ B gav k _ Q.7668 X 0-036 X 0 0254 X 1.5 
ii„H m 9 ' 1.26 X 10"* X 71.99 X 10 3 


The emf equation is: 


(5-1) 

(5-2) 


(5-3) 
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E m = B,„ W^'t^nDi 


(5-4) 


where IV, is the number of turns on stator. Thus, from (5-4) with IV, = 21 turns. 


A- 


E m 


240 


*B g „ W 1 u l e4 nkt X X 0.7668 X 21 X 628 X 1 X 0.0254 
= 0.2963 m = 11.62 in 


(5-5) 


5.4. Summary of Design Data 

A = 11 in = 0.2805 in 
A m = 0.46 in 
g = 0.036 in 
t = 1.2 :n 

T m = 1 *n 

9m = T-T m = 0.2 m 


5.5. Further Calculations 


K ■■ 


T 

9m 


(A — 2j — A m )A m + — h~ 


+ 2(A - 2g - A m ) 


+ (A - 9 ~ A m )ln 


A m + -9 




Leakage flux of permanent magnet is: 

4>„ = — KH m = - KH C (1 - ^-) 

B, 

The flux linking with the stator coil is: 
4 > — *(A — - A m )r m B m - 

which is related to the induced voltage by: 

£ = V2 7T W, / m 0 

With the current density chosen as J (9 
given by: 


(5-7) 

(5-8) 

(5-9) 

= 5 A/mm 2 , the wire cross-section area is 


A co =-j- = = 41.66 mm 2 = 0.0641 in 2 

Jco 5 
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To obtain this area of cross-section, 4 No. 7AWG wires, each having a diameter 
of 0.1443 in, may be used in parallel. Thus: 


A eo = 4 X — (0.1443 X 25.5) 2 = 42.54 mm 2 
4 

Let the slot fill-factor k M = 0.5. Then slot area: 

A = , 21X4234 _ mm , _ 2J5 




0.5 


(5-10) 


Let b t = 0.5 in, = 0.45 in, Ai =0.2 in, b p = 0.7 in, A 2 — 0.8 in, and b, 2 = 3.44 in. 
Then: 


b, 2 h 2 — 3.44 X 0.8 = 2.753 in 2 = A (5-11) 

which is adequate slot area. 

By the computer calculation for above design, the following no-load data are 
obtained: 


E = 169.505 V 
<j> = 0.01817 W b 
= 1-1446 
B m = 0.9724 T 


= -70.252 M/m 

Average diameter of the winding: 


D„ = A + 2Aj + A 2 = 11 + 0.4 + 0.8 = 12.2 tn 


Stator winding resistance: 


A. = P 


= 2.1 X 10 -8 — - X ir —" 12 ' 2 X 0 0254 = 0.0103 Q 
A c<> 42.54 X 10"* 


(5-12) 

(5-13) 


Stator leakage inductance: 


L„ = ft'irWf 


(Pi + 2Aj + 2hn)h 2 2 A A) 

3 A, 2 b, i + b, 2 


Af 2A, 2 
2A,2 A,i + b, 2 


= 1.26 X 10“* X ff X 21 2 


13 X 0.8 
3 X 3.44 


2 X 11 X 0-2 
2.04 + 3.44 


0.8 2 2 X 0.2 2 

2 X 3.44 + 2.04 + 3.44 


X 0.0254 


= 0.0771 mH 


(5-14) 


For 7= 208.33 A and 6 = 45°, equivalent magnetizing inductance L mt = 0.832 
mH, as calculated by computer. This inductance accounts for armature reaction as 
well as the effect of change of the magnet leakage flux. Hence, the total synchronous 
reactance of the linear alternator becomes: 

X, = w r (L mt + L„) = 628(0.832 + 0.0771) X 10" 3 


= 0.5709 Q 


(5-15) 
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5.6. Maximum Power Condition 

Suppose that the load of the linear alternator is purely resistive and armature 
resistance is negligible, the phasor diagram of the linear alternator is shown in Fig. 
5-1. 


Output power: S = VT 

(5-16) 

From Fig. 5-1, 


/ X, =* j? sin 0 

(5-17) 

or 


/ => sin0 

-X* 

(5-18) 

And 


V — E cos 0 

(5-19) 

Thus, 


.E 2 E 2 

5 = — sin 0 cos 0 = —— sin20 

2X t 

(5-20) 

which is maximum when 20 = 90° or 0 = 45'. Under this condition, 


r E 169.505 . 

vS*" Vi X 0.5709 - 209 - 95 * 

(5-21) 


Rated Load: R L = X, - R a « 0.5709 - 0.0101 = 0.5608 D 
Output Voltage: V = IR L = 209.95 X 0.5608 = 117.74 V 
Output kVA = 117.74 X 209.95 X 10 -3 = 24.72 kVA 


5.7. Weight-to-Power Ratio 

( Vol )fUr = 7 [(A + 2 A t + 2 h 2 + 2 b b f - (£>,. + 2 A, + 2 A 2 ) 2 ]( 2 & t + 6 ^) 
+ f [(A + 2A, + 2A 2 ) 2 - (A + 2A,) 2 ] 2b t 


2 7T 


[(A + 2A) - Df\ 


2 (bt + r) 


= J (13.9 2 - 13 2 ) X 3.54 + f (13 2 - 11.4 s ) X 1-0 + - (11.4 2 - ll 2 ) X 0.85 


= 109.936 in 3 = 0.001823 m 3 

(™Wu. = 4 [(A - 2?) 2 - (A - 2g - 2h m f\ j- r m 
= (10.94 2 - 10.02 2 ) XJTXl 


(5-25) 


= 60.58 in 3 = 0.001005 m 3 


(5-26) 
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(™W, = 7 [(A - 2? - 2A m ) 2 - (£>,• -2j- 2A m - 2A p ) 2 ](6, 2 + 2r) 

= 7 (10.02 2 - 8.62) 2 (3.44 + 2.4) 

= 119.695 m 3 = 0.001985 m 3 (5-27) 

( V °l)" PP 'r = » A, W, 

= Jr X 12.2 X 0.06542 X 21 = 52.655 m* = 0.000873 m 3 (5-28) 

Assuming a mean mass density of 7.8 x 10 3 kg/m 3 , the total weight of the 
linear alternator becomes: 


W = (0.001823 + 0.001005 + 0.001985 + 0.000873) X 7.8 X 10 3 
= 44.3508 kg = 97.572 lb. 
weight _ 44.3808 


power 


24.72 


1.794kg AW = 3.947 lb AW 


(5-29) 

(5-30) 


5.8. Losses and Electrical Efficiency 

Copper loss: 

pR t = 209.95 2 X 0.0101 = 445.198 W 
Core loss: 

First we determine the weight of iron: 

= 2 X 7 [(A + A,) 2 - A 2 ] J ( b t + t)x 7.8 X 0.254 3 
= 1.5473 kg = 3.404 lb 

= 7 [(A + 2 A, + 2 A 2 ) 2 - (D ( + 2h x f^b t X 7.8 X 0.254 3 


= 3.9656 kg = 8.724 lb 


^ (A + 2A[ + 2A 2 + 2 6$) 2 — 


— (A + 2Aj + 2A 2 ) 2 


(A ,2 + 26 t ) X 7.8 X 0.254 3 


(5-31) 


(5-32) 


(5-33) 


= 8.7056 kg = 19.153 /A 


(5-34) 
(5-35) 

From the flux density of every port of iron and core loss curve shown in Fig. 5- 
2, the core loss coefficients are found as follows: 


W p = 1.985 X 7.8 = 15.483 kg = 34.063 lb 


0.304 


100 


60 


|1.3 


, A,, = 0.587 


100 


M.3 


60 


i Ac* = 0-609 


t 'll.3 
100 


60 


and A = 0.957 


100 


V.3 


60 
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Figure 5-2 Core loss and exciting VA for M-19 steel at 60 Hz. 
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The corresponding losses are: 


t V 


= 3.404 X 0.304 • 

100 

60 

W K ,C K , = 8-724 x 0.587 

100 

60 

W Kh C K * = 19.153 X 0.609 

i 100 
( 60 

W,C P = 34.063 X 0.957 ■ 

100 I 1 
60 


|1.3 


2.0103 W 


1I.3 


= 9.949 W 


V-3 


22.66 W 


1.3 


63.329 W 


{lo»s) Mal „„ = (loss)„ + (lo$s) K , + (loss) Kh + (loss), 

= 97.948 W 

copper loss + core loss = 445.198 + 97.948 = 543.146 W 


Electrical efficiency =- output - 

output + losses 


24.72 


24.72 + 0.543 


= 97.85% 


(5-36) 

(5-37) 

(5-38) 

(5-39) 

(5-40) 

(5-41) 

(5-42) 
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CHAPTER 6 

OPTIMAL DESIGN OF A PM LINEAR ALTERNATOR 

6.1. Concept of Optimal Design 

For a specified rating (kVA) and output voltage (or current), a PM linear alter¬ 
nator can have numerous possible shapes and dimensions. The objective of optimal 
design is to design the machine which has the best pre-determined characteristics and 
satisfies specified constraints such as output voltage and alternator rating. In our 
design, we will attempt to achieve a minimum power-to-weight ratio and a max¬ 
imum electrical efficiency by applying a mathematical optimization technique. 

In an optimal design, we choose design variables and an objective function, 
which is a function of design variables and describes the characteristics to be optim¬ 
ized. We then find the set of values of the design variables, which correspond to the 
maximum or minimum value of the objective function with some equality or inequal¬ 
ity constraints, which are the required conditions that the design must satisfy. 
Mathematically, the optimization problem may be stated as follows: 

For the design variables x It x 2 ,...,x n = X, find 


min F = f(x v x 2 ,...,x„) = /(X) 

(6-1) 

subject to: 


9j{X) <0 j = 1,2,...,/ 

(6-2) 

A t (X) = 0 k = 1,2 

(6-3) 


6.2. Penalty Function Methods 

The penalty function methods convert a constrained optimization problem into 
an unconstrainted optimization problem. These methods are algorithms which gen¬ 
erate a sequence of points from X< 0) to X< r) , where X< 0) is the initial point, X*' 1 the gen¬ 
eric point and X^ the limit point, and the best estimate of the optimal point X* is 
produced by the algorithm. The points X^, f-l,2,...,7 are stationary points of the 
penalty function. 

Consider the penalty function: 

P(X,R) = f(X) + J?[fl,*(XU(X)] (M 

where R is a set of penalty parameters, and 0, the penalty term. Then P{X,R) is a 
function of R and the constraint functions. Penalty function methods are classified 
according to the procedure employed for handling inequality constraints. These 
methods are referred to as interior or exterior point method according as the 
sequence X^ contains feasible or infeasible points, respectively. If the sequence of 
stationary point contains both feasible and infeasible points, the method is said to be 
mixed. 

The interior penalty function technique arguments the objective function with a 
penalty term which is small at points away from the constraints in the feasible 
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region, but which "blow up" as the constraints are approached. 

The penalty function is of the form: 

p(F,R t) = ,m - h ^ | *. m («) 

with R t > R 2 > • • • > R t > Rt +1 > • • • >0 
and lim R t = 0 

t—*O0 

Therefore, when t -* oo, —^—► oo, and it forces A*^'*] 2 , k=l,2,...,K, approach zero, 

when minimum of P{F,R t ) is approached. Inequality constraints are satisfied since 
t=l,2,...,T, are feasible points of inequality constraints. When t—*x>, P{F,R t ) has 
the same external as F =* f{X) has. 

The exterior penalty function technique modifies the objective function by 
adding a penalty whenever a constraint is violated. The penalty function takes the 
form: 

P(F,R t ) = /[*«] + Rt t < f/W > 2 + E (8-6) 

y-1 kM 

with 0 < R x < R 2 ... < R t < R t + i < •— 

lim R t == + oo 

MX) 


and 


<*OT> 


o if ?< [^‘)] < o 

\g { m if > o 


With t increasing, the minimum of P(F,R t ) will conjugate to the minimum of 
F = f{X) and the constraints of g,{X] will be satisfied for Also the A t p^] 2 

will be force to approach zero when a minimum of P(F,R t ) is approached for 
k—1,2,...,K. 


6.3. Optimal Design of PM Linear Alternator 

For the PM linear alternator shown in Fig. 2-1, we choose: r m = 1.0 in, A, — 0.2 
in. Let: = g\ X 2 = A m ; X s = D { ; X i = r; X 8 = X t = 6 P ; X 7 = A 2 ; X 8 = i, 2 ; X 9 = W l 

(number of turns of the coil); and X x0 =» b t are the design variables, which change dur¬ 
ing the optimization procedure. Other dimensions can be obtained from r m , A! and 
the design variables. The cross-section of the wire used for the coil is given by: 

A„ - (8-7) 

where k m is the slot fill factor = 0.5 (assumed). 

The objective function in this design is the sum of the ratio of weight to output 
power (Ib/kVA) and the ratio of iron and copper losses to output power, i.e.: 

p _ Ei 

S + 10005 


( 6 - 8 ) 





where W - weight in lb; S = kVA rating; Z { — iron loss in watts; and Z c = copper loss 
in watts. 

The goal of the optimal design is to find a set of design variables which make F 
a minimum and satisfy the following constraints: 

output power S’ = 25 kVA 
output voltage V < 500 V 
current density J c < 5A/mm 2 
pole shoe length r < 1.2 in 
airgap g > 0.03 in 

The optimization subroutine package OPTLIB was used to obtain the optimal 
design, and the method used is the exterior penalty function. So, the optimal design 
problem reduces to obtaining the design variables which correspond to the minimum 
value of an unconstrained function P(F,R), which is a function of F and the above 
constraints. The algorithm of optimal design is as follows: 

(i) Choose initial values of design variables. 

(ii) Calculate no-load electromagnetic field, no-load induced voltage E 0 , coil 
resistance R t and leakage inductance L,„. 

(in) Calculate electromagnetic field corresponding to 9 * 45* and purely resistive 
load. Find the induced voltage drop due to armature reaction and the work¬ 
ing point variation jf the permanent magnet. Then find the equivalent 
mutual inductance L me , output voltage V, stator current I u output power 5, 
weight W and the losses Z { and Z e . 

(iv) Check if P{F,R) has a minimum value and satisfies accuracy 

yes: Go to (vi) 
no: Go to (v) 

(v) Find new search direction and new design variables, go to step 2. 

(vi) Stop. 

The flow chart of optimal design is shown in Fig. 5-1. 

The initial values of the design variables were chosen as follows: 
aq * g = 0.03 in 
x 2 — h m = 0.4 in 
x 3 = D { = 16.0 in 
x i — t — 1.1 in 
= & 6 = 0.8 in 
= 6 p = 0.8 in 
x 7 — hn = 2.0 in 
Xg “ b t 2 “ 3.5 m 
x 3 = Wj = 21 turns 
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Figure 6-1 Flow chart of optimal design. 
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The optimal values of the design variables are: 
x x = g = 0.031 in 
x 2 = h m *■ 0.44 in 
x 3 = D- = 11.35 m 
x 4 = r = 1.2 in 
x 5 — h ** 0.34 m 
x fl = b p — 0.50 in 
x 7 = A 2 = 0.91 in 
Xg = 6,2 “ 3.63 in 
x 9 = Wj =ss 85 tarns 
xjo = — 0.38 in 


From the design the following characteristics of the linear alternator 
obtained: 

No-load induced voltage, E„ = 571.58 Volt 

Output voltage, V = 395.26 Volt 

Stator current, 7 , = 62.74 A 

Stator resistance, R t — 0.142 f2 

Stator leakage inductance, L, a = 1.93 mH 

Stator mutual inductance, L m = 8.33 mH 

Cross-section of stator winding A t0 = 12.715 mm 2 

Weight/power ratio, = 3.72 Ib/kVA 

Losses/output power Z/1000S = 0.0233 
Output power, S = 24.8 kVA. 

Field of without load: 


B m - 

B u , = 

= 

Bpek = 

Field for 7, = 62.74 A and Q = 45° 

B m = 0.82138 T 

B„ = 0.48215 T 

B , c , = 0.90188 T 

B Kk = 0.94491 T 

Bp,. = 0.73906 T 

= 1.51003 T 


H m = 

-129306.30000 A/m 

= 

66.29997 A/m 

h k , = 

2588.85600 A/m 

H, eh - 

3974.03300 A/m 

77p„ 

83.88339 A/m 

= 

15728.46000 A/m 


= 

-178987.50000 A/m 


38.75073 A/m 

- 

77.73460 A/m 

- 

83.11403 A/m 

= 

59.51506 A/m 

flpc* - 

2477.14200 A/m 


0.89039 T 
0.81040 T 
1.51587 T 
1.58820 T 
0.95107 T 
1.94321 T 


are 
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CHAPTER 7 
CONCLUSIONS 

In the preceding chapters certain important aspects of analysis, design and 
optimal design of single-phase single-slot tubular PM linear alternators have been 
investigated. 

Methods developed in Chapters 2 and 3 are used to find the mutual magnetic 
field in every section of the linear alternator and the leakage magnetic field, with 
load or without load. Furthermore, the induced voltage may be obtained from the 
magnetic field, if the stator winding turns and frequency of the alternator are known. 

The field of every section of the alternator, found by the methods developed in 
Chapter 2 and 3, is the average value of the field. However, it is good enough for 
finding the induced voltage. The advantage of these methods is that they make the 
calculation of the alternator’s field much more easy and simple than by numerical 
method. The simplicity is significant for the optimal design of the alternators. 

The methods and formulas presented in Chapter 4 are available for finding the 
parameters of the alternators, which are necessary for transient state and steady 
state analysis of the alternators. 

The procedure discussed in Chapter 5 is the essential design guideline of tubular 
PM linear alternators. Combining the design guideline in Chapter 5 and mathemati¬ 
cal optimization method, the technique presented in Chapter 6 ensures that the 
designed alternator has best specialized characteristics, which is defined by the 
optimal function, as shown in the design example. 

The methods and formulas presented in this report and in reference [4] and [5] 
form the essential knowledge for analysis, design and operation of PM linear alterna¬ 
tors. 
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C 


HO-O.5*HC 
10 CONTINUE 

BM—BR*(1.0-H0/HC) 

QL—AK*H0 ^ 

Q-PI*(DI-2.0*GG-HSM)*TM*BM-QL 

BSP-Q/SSP 

BSCV-Q/SSCV 

BSCH-Q/SSCH 

BPCV-(Q+QL)/SPCV 

BPCH-(Q+QL)/SPCH 

CALL BH(BSP,HSP) 

CALL BH(BSCV,HSCV) 

CALL BH(BSCH,HSCH) 

CALL BH(BPCV,HPCV) 

CALL BH(BPCH,HPCH) 

BGAV-Q/SG1 

HGAV-BGAV/U0 

FM-HSCH*B+HSCV*(2.0*H2+BB)+HSP*2.0*H1+HPCV*BP+HPCH*B+HGAV*2.0*GG 
HM—FM/2.0/HSM 

IF (ABS(HM-H0).LT.1.0) GO TO 50 
H0-0.5*(HM+H0) 

GO TO 10 

50 HM—0.5*(HM+H0) 

BM-BR-BR/HC*HM 

AKM-(Q-AK*HM)/Q 

E0-SQRT(2.0)*PI*W1*FME*Q 

AS-H2*BS2 

AW-0.5+AS/W1 

RA-2.1E-8*W1*PI*(DI+2.0*H1+H2) /AW 
BS-0.5*(BS1+BS2) 

ALL—PI*W1**2*U0*((DI+2.0*(H1+H2))*H2/3.0/BS2-H2**2/2.0/BS2+ 

* 2.0*(DI+H1)*H1/BS) 

AIN-1.0E6*AW*5.0 
AI0-AIN 

D-l.0-AK*HC/BR/SG3M 

D1-BR*SG3M*GG/SG1 

D2-SG3M*GG/SG1/SG3 

A—Dl/ ( (GG+HSM) /SG2+D2-GG/SG11) 

C-(-U0*HSM+Dl/HC*D)/((GG+HSM)/SG2+D2-GG/SG11) 

105 AIF-AI0/SQRT(2.0) 

HO-O.5*HC 
110 CONTINUE 
Q2-A+C*H0 

Ql-((HC-D+HO)*BR/HC+Q2/SG3)*SG3M 

Q-Q1-Q2 

BSP-Q/SSP 

BSCV-Q/SSCV 

BSCH-Q/SSCH 

BPCV-(Q1-AK*H0-Q2)/SPCV 
BPCH-(Q1-AK*H0-Q2)/SPCH 
CALL BH(BSP,HSP) 

CALL BH(BSCV,HSCV) 

CALL BH(BSCH,HSCH) 

CALL BH(BPCV,HPCV) 

CALL BH(BPCH,HPCH) 

FM—HSCH*B+HSCV*(2.0*H2+BB)+HSP*2.0*H1+HPCV*BP+HPCH*B 
HM-((-Q2/(Q1-Q2)*FM+W1*AIF)*SG2*U0/2.0/(GG+HSM)-A)/C 
IF (ABS(HM-HO).LT.1.0) GO TO 150 
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i 

i 


C 


c 


c 


H0-0.5*(HM+H0) 

GO TO 110 
150 HM-0.5*(HM+H0) 

BM=BR-BR/HC*HM 

Q2-A+C*HM 

Ql-((HC-D*HM)*BR/HC+Q2/SG3)*SG3M 
Q-Q1-Q2 

AKM-(Q1-AK*HM)/Ql 

E-SQRT(2.0)*PI*W1*FME*Q 

AL»(E0-E)/AIF/628.0 

AI-E0/SQRT(2.0)/(AL+ALL)/62 8.0 

IF (ABS((AI-AI0)/AI0).LT.0.03) GO TO 152 

AI0-0.5MAI+AI0) 

GO TO 105 
152 CONTINUE 


Dl-DI+2.0*(H1+H2+BB) 
D2—DI+2.0*(H1+H2) 
D3-DI+2.0*H1 
D4-DI-2.0*GG 


D5-DI-2.0*(GG+HSM) 

D6-DI-2.0*(GG+HSM+BP) 

Vl-PI/4.0* (Dl**2-D2**2) * (BS2+2.1)*BT) 
V2-PI/4.0*(D2**2-D3**2)*2.0*BT 


V3-PI/4.0* (D3**2-DI**2)*0.5*(BT+T)*2 0 
VS-V1+V2+V3 


VM-PI*(D4**2-D5**2)*TM 

VP-PI/4.0*(D5**2-D6**2)*(BS2+2.0*T) 

VC-PI*(DI+2.0*H1+H2)*AW*W1 

W-7.8*(VS+VM+VP+VC)*1000.0*2.2 

V-E0/SQRT(2.0)-AI*RA 

S-V*AI 


ZC-AI**2*RA 
CALL CL(BSCH,FSCH) 

CALL CL(BSCV,FSCV) 

CALL CL(BSP,FSP) 

CALL CL(BPCH,FPCH) 

FRA-(100.0/60.0)**1.3 
ZSCH-2200.0*V1*FSCH*FRA 
ZSCV-2200.0*V2*FSCV*FRA 
ZSP-2200.0*V3*FSP*FRA 
ZP-2200.0*VP*FPCH*FRA 
ZR-ZSCH+ZSCV+ZSP+ZP 
- Z-ZC+ZR 

F-(1000.0*W+Z)/S*10.0 
G(6)-(25000.0-S)/25000.0 
G (3)-(AI-AIN)/10.0 
G(2)—(0.3-X(1))*100. 
G(1)-(X(4)-0.12)*10. 

G(4) — (1.6 *RR-BS1)*10 

G(5)-V-500.0 

UF-F 

RETURN 

END 


SUBROUTINE GRADU(X, S, SUM) 

RETURN 

END 


SUBROUTINE CL(B,F) 
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339 
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342 

343 

344 

345 
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347 


F-(B-0.4)*1.9/1.2+0.1 
RETURN 
END 
C 
C 

SUBROUTINE BH(B,H) 

DIMENSION BC(12),HC(12) 

DATA BC/0•0,0.05,0.2,0.4,0.6,0.8,1.0,1.2,1.4,1.6,1.75,1.8/ 
DATA HC/0.0,13.0,22.0,33.0,47.0,65.0,90.0,130.0,370.0, 

* 4200.0,8000.0,10000.0/ 

DO 130 1-2,12 

IF (B.LT.BC(I).OR.B.GT.BC(12)) GO TO 110 
GO TO 130 

110 IF (B.GT.BC(12)) GO TO 120 
L-I-l 
GO TO 140 
120 L—11 

GO TO 140 
130 CONTINUE 

140 H-HC(L)+(B-BC(L))*(HC(L+l)-HC(L))/(BC(L+l)-BC(L)) 

RETURN 

END 

C 

..INC OPTLIB FORTRAN 
C 
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